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Abstract. - We demonstrate the appearance of a purely entropic stochastic resonance (ESR) oc- 
curring in a geometrically confined system, where the irregular boundaries cause entropic barriers. 
The interplay between a periodic input signal, a constant bias and intrinsic thermal noise leads to 
a resonant ESR-phenomenon in which feeble signals become amplified. This new phenomenon is 
characterized by the presence of two peaks in the spectral amplification at corresponding optimal 
values of the noise strength. The main peak is associated with the manifest stochastic resonance 
synchronization mechanism involving the inter-well noise- activated dynamics while a second peak 
relates to a regime of optimal sensitivity for intra-well dynamics. The nature of ESR, occur- 
ring when the origin of the barrier is entropic rather than energetic, offers new perspectives for 
novel investigations and potential applications. ESR by itself presents yet another case where one 
constructively can harvest noise in driven nonequilibrium systems. 
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, Stochastic resonance (SR) is an intriguing phenomenon 
occurring in systems pertaining to the wide class of pe- 
r,iodically modulated noisy systems which includes many 
different theoretical and experimental situations [1-10]. 
SR relates to a remarkable idea that changed our com- 
mon perception of noise [4]: in particular, ambient noise 
ipay play a constructive role in amplifying feeble signals 
or may facilitate noisy transport. The phenomenon has 
been well studied mainly in systems having an intrinsic 
energetic potential whose origin is the presence of inter- 
4ctions. However, the prevalent role of interactions in the 
dynamics of a system is by no means general. It is indeed 
the free energy what controls the dynamics and it could 
happen that the entropic contribution plays a leading role. 
One example is the case of ion channels, where the role of 
entropic contributions for the phenomenon of SR has also 
been addressed [11-14]. 

Here, our focus is on stochastic resonance phenomena in 
confined systems, where stylized, purely geometrical con- 
strictions may lead to a dominant entropic potential with a 
strong impact on the transport characteristics [15-18] that 
may exhibit a SR behavior in some situations [19-21]. In 
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this Letter we discuss a purely entropic resonant behavior 
in a geometrically confined system in the presence of a lon- 
gitudinal, constant bias and an oscillating force. This new 
resonant phenomenon is distinctly different and goes be- 
yond the standard SR picture in several aspects. First, its 
origin is strictly entropic and is solely associated with ge- 
ometric unevenness and confinement; in particular, there 
is no energetic barrier in the system. Second, the situa- 
tion is characterized by the presence of two peaks in the 
amplification factor signaling two different optimal values 
of the noise strength. Finally, the enhancement of the 
amplification associated to the second peak is not due to 
noise activation, but rather is due to an optimal sensitiv- 
ity of the intra-well dynamics to noise; it emerges when 
the entropic barrier disappears. 



To illustrate this phenomenon, we consider the dynam- 
ics of a Brownian particle in the two-dimensional struc- 
ture depicted in fig. 1. The time evolution of this particle, 
occurring in a constrained geometry subjected to a sinu- 
soidal oscillating force F{t) and a constant bias Fb acting 
along the longitudinal direction of the structure, can be 
described by means of the Langevin equation, written in 
the overdamped limit as 
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Fig. 1: Schematic illustration of the two-dimensional structure 
confining the motion of the Brownian particles. The symmet- 
ric structure is defined by a quartic double well function, cf. 
eq. (2), involving the geometrical parameters L^, Ly and b. 
Brownian particles are driven by a sinusoidal force F{t) and a 
constant bias _Fb acting along the x— direction. 
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dr 

dt 



F{t) + ^jkBT^it), 



(1) 



where r denotes the position of the particle, 7 is the 
friction coefficient, Ex is the unit vector along the ax- 
ial direction x, and ^(t) is a Gaussian white noise with 
zero mean which obeys the fluctuation-dissipation rela- 
tion (CjW6(^')> = 2S^jS{t-t') for i,j = x,y. The 
explicit form of the periodic input signal is given by 
F{t) = Asm{nt), where A is the amplitude and fl is 
the frequency of the sinusoidal signal. In the presence 
of confining boundaries, this equation has to be solved by 
imposing no-flow boundary conditions at the walls of the 
structure. 

The shape of the considered 2D structure, fig. 1, which 
is mirror symmetric about its x— and y— axis, is defined 
by 



WI^X) = Ly 



X 



-2i„ 



-2 =-Mx)^ (2) 



where wi and uiu correspond to the lower and upper 
boundary functions, respectively, corresponds to the 
distance between the location of the bottleneck and that of 
maximal width, Ly quantifies the narrowing of the bound- 
ary functions and b is the remaining width at the bottle- 
neck, cf. fig. 1. The local width of the structure is given 
by 2w{x) = Wu{x) — wi{x). This particular choice of the 
geometry is intended to resemble the classical setup for 
SR, namely a double well potential, where two basins are 
separated by a barrier. However, we shall show that the 
distinct nature of this confined system, dominated by a 
purely entropic rather than energetic potential landscape, 
gives rise to new phenomena. 

For the sake of simplicity, we use reduced, dimensionless 
units. In particular, we scale lengths by the characteristic 
length ij;, time by r = ^iL^^ /k^T^, which is the char- 
acteristic diffusion time at an arbitrary, although typical 
reference temperature Tr, force by Fr = k^T^/Lx, and 
frequency of the sinusoidal driving by 1 /r. In dimension- 
less form the Langevin equation (1) and the boundary 



functions (2) read: 



dr 
di 



= F^,ex-F{t) ex + VD^{t), 



wi{x) = — Wu(.t) = —w{x) = ex^ — 2ex^ — b/2 . 



(3) 
(4) 



where D = T/T^ is the dimensionless temperature and 
e = Ly/L^ is the aspect ratio of the structure. 

Reduction of dimensionality. In the absence of 
a time-dependent applied force, i.e., F{t) ~ 0, and tilting 
force, i.e., i^b = it has been elaborated previously [15-20, 
22] that the 2D Fokker-Planck dynamics corresponding to 
the Langevin equation (3) can be reduced to an effective 
ID Fokker-Planck equation, reading in dimensionless form 



dP{x,t) 
di 



d f^ dPjx^t) 
dx 1 dx 



V'{x,D)P{x,t) 



(5) 



where the prime refers to the derivative with respect to x 
and the effective potential [15,20] 



V{x,D) = -D \n[2w{x)] . 



(6) 



This equation describes the motion of a Brownian par- 
ticle in a bistable potential V{x^ D) of purely entropic 
nature. It is important to emphasize that this bistable 
potential was not present in the 2D Langevin dynamics, 
but arises due to the geometric restrictions associated to 
confinement. In general, the diffusion coefficient will de- 
pend on the coordinate x as well [16,17,22-24], but in the 
case discussed here this correction does not significantly 
improve the accuracy of the reduced kinetics. Therefore 
wc do not further consider this correction. Notably, for 
vanishing width 2w{x) at the two opposite corners of the 
structure in fig. 1 this entropic potential approaches in- 
finity, thus intrinsically accounting for a natural reflecting 
boundary. 

In the presence of a constant bias along the x— direction 
of the channel the free energy becomes [15-17] 



V{x, D) = -Fb x-D ln[2 w{x)] . 



(7) 



The behavior of this tilted potential is depicted in fig. 2. 
In contrast to the purely energetic potentials used for 
classical SR [1,25], the barrier height and the shape of 
the potential change with noise strength D. At small D 
the strength of the entropic contribution is almost negligi- 
ble and the potential can be approximated by V{x, D) « 
— i^b X. The corresponding one-dimensional Fokker-Planck 
equation in the absence of periodic force reads 



dP{x,t) 

m 



d_ 

dx 



D 



dP{x,t) 
dx 



F^P{x,t) 



(8) 



However, in general, the potential function V{x,D), 
eq. (7), describes a tilted bistable potential, where the 
height of the potential barrier depends on the noise 
strength D. The location of the minima also depends on 
the noise level, tending to a limiting value of .t = ±1 for 
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Fig. 2: (Color online). The effective one-dimensional potential 
V{x,D), eq. (7), at a constant bias Fh and for increasing noise 
strength D. The aspect ratio e and the bottleneck width b of 
the structure are set to 1/4 and 0.02, respectively. The grey 
shaded regions on the left and right resemble the impenetrable 
boundaries. The potential function diverges at and xn. 



large D. More importantly, the barrier disappears and the 
potential does not exhibit any minima when D is smaller 
than a critical value Dc, which depends on the geometrical 
structure and on the tilt. For the shape given by eq. (2) 
with e = 1/4 and b — 0.02, the potential has an inflection 
point (i.e., the barrier disappears) when Dc ~ 0.144i^b- 
We will discuss below the relevance of this special situa- 
tion in the context of SR. 

Spectral amplification. — The response of the sys- 
tem to the weak periodic input signal, F{t) — Asui{Qt), 
is also a periodic function of time. The mean value of the 
position x{t) in the asymptotic time limit (i.e., after the 
memory of the initial conditions is completely lost) admits 
the Fourier series representation [26, 27] 



{x{t))as 



E 



I at 



(9) 



with the complex-valued amplitudes M„ which depend 
nonlinearly on the driving frequency, driving amplitude 
and on the noise strength. 

In order to study the occurrence of SR we analyzed the 
response of the system to the applied sinusoidal signal 
in terms of the spectral amplification rj, which is defined 
as [26,27] 



V = 



2\Mi\ 
A 



(10) 



Tvifo-State approximation. — If the potential bar- 
rier is sufficiently high and the two basins of attraction are 
well separated, the intra-well motion may be neglected, 
and a simplified two-state description can be used to get 
useful insights into the full dynamics [28,29]. Within 
the two-state approximation the transition rate from one 
well to the other can be determined in terms of the mean 
first passage time (MFPT) [25] to reach a potential min- 
imum starting out from the other minimum of the tilted 



bistable potential V{x,D). Accordingly, the forward (fc_|_) 
and backward (fc_) rates are given by 



T±{x^ x± ) 



(11) 



where x- and x+ indicate the noise(£')-dependent loca- 
tion of the left and right minima, respectively; T+(x_ — > 
a;+ ) is the mean first passage time for reaching 2;+ starting 
out at x-] and, vice- versa, T-{x+ — > a;_) is the mean first 
passage time in the backward direction. More explicitly, 
the forward and backward mean first passage times are 
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wix) 



dyw{y)e''''y^° , (12) 



and 



-Fbx/D px 



dx ■ 



w{x) 



dyw{y)e^^y/° , (13) 



where xl and are, respectively, the left and right lim- 
iting values at which the boundary function vanishes, see 
fig. 2. Note that the bias i^b must be small in order to have 
two well separated basins of attraction (see discussion of 
the potential) . Within the two-state description the mean 
position is given by 



(14) 



A sinusoidal signal leads to a modulation of the barrier 
height of the double well potential and consequently to a 
modulation of the rates k± [28]. One obtains the spectral 
amplification in the lowest oder of A/ D, reading 



V = 



1 



{x+ — X-)'^ fc_ 



1 



(fc+ + A:_)2 + J72 



(15) 



Note that for the symmetric case, i.e., = = k and 
x± = ±1, the spectral amplification given in eq. (15) re- 
duces to the well known expression. 



4fc2 



V 



1 



D2 4fc2-|-f72 



(16) 



which in the linear response regime is by construction in- 
dependent of the amplitude strength [1,26,27,29]. 

ID modeling. — The two-state description of the 
linear response fails for either large driving amplitudes, 
small driving frequencies [1,30] or very weak noise [31]. In 
such circumstances, the ID Fokker-Planck equation has to 
be considered and the nonlinear response analyzed [1,27, 
30,31]. In the presence of a sinusoidal signal F{t) the ID 
Fokker-Planck reads 



dP{x,t) 
Ft 



+ F(t)]P\. (17) 



We have integrated numerically the above equation 
by spatial discretization, using a Chebyshev collocation 
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Fig. 3: (Color online). The time averaged mean position Mq 
as a function of the noise strength D for different tilting force 
values. The solid lines correspond to the numerical integration 
of the ID Fokker-Planck equation, eq. (17), and the dashed 
lines correspond to the two-state approximation, eq. (14), i.e.. 
Mo = qo. In the potential function we set the aspect ratio 
e = 1/4 and the bottleneck width b = 0.02. The dotted lines 
represent the behavior of Mo at small D obtained by integra- 
tion of eq. (8). For the symmetric case (Fb = 0) all these lines 
collapse to the black solid line, i.e. Mo = 0. 



Fig. 4: (Color online). The behavior of the spectral amplifica- 
tion factor as a function of the noise strength D for the same 
parameters as in fig. 3. The solid lines correspond to the nu- 
merical integration of the ID Fokker-Planck equation, eq. (17), 
whereas the dashed lines correspond to the two-state approx- 
imation, eq. (15). The behavior of rj at small D obtained by 
integration of eq. (8) with a periodic forcing and at the respec- 
tive tilting force values is represented with dotted lines. The 
arrows indicate the critical value Dc that corresponds to the 
vanishing of the barrier. 



method, and employing the method of lines to reduce the 
kinetic equation to a system of ordinary differential equa- 
tions, which was solved using a backward differentiation 
formula method. Thereby wc obtain the time-dependent 
probabihty distribution P{x,t) and the time-dependent 
mean value of the position, {x{t)) = J x P{x, t) dx. In 
the long-time limit this mean value approaches the peri- 
odicity of the input signal with angular frequency , see 
in Refs. [26,27]. After a Fourier-expansion of {x{t)) one 
finds the time averaged mean position Mo = Mo{A, D) 
and amplitude Mi ~ Mi{A, D) of the first harmonic of 
the response. 

The mean position Mq as a function of the noise 
strength D is depicted in fig. 3 for two different tilt- 
ing force values. For comparison, the approximated mean 
position qo obtained within the two-state modeling is also 
plotted. At very small D there is no barrier and the par- 
ticle oscillates just in the vicinity of the boundary, either 
at xr or Xh depending on the direction of the tilt. As the 
noise strength increases the barrier height increases, see 
fig. 2, and the inter- well (from one well to the other) dy- 
namics become more dominant. Remarkably, the behavior 
of the mean position can be qualitatively captured within 
the two-state approximation. 

For the ID modeling the spectral amplification for 
the fundamental oscillation is computed accordingly, cf. 
eq. (10). The comparison between the results of the ID 
modeling and the two-state approximation for the param- 
eter ry as a function of the noise strength D and for two 
different tilting force values is depicted in fig. 4. Within 
the ID modeling 77 exhibits a double-peak behavior for a 
finite tilt. The appearance of the main peak at higher D, 
which is due to the synchronization of the periodic signal 



with the activated inter-well dynamics, can be nicely cap- 
tured within the two-state model. However, there is a sec- 
ond peak at small values of noise that cannot be described 
within the two-state model. As it is discussed before, at 
small noise strengths there is no barrier and the particle 
oscillates in the vicinity of the boundary. As we increase 
the noise one can observe a steep rise in 77 with D which is 
a consequence of noise helping the particle to climb higher 
the potential hill (see fig. 4). The second peak at small 
D is attributed to the intra-well dynamics, and could also 
be observed for mono-stable energetic potentials [5,32,33]. 
The steepness of increase of 77 depends on the strength of 
the tilting force. However, beyond D^, the barrier appears 
and gets higher as wc increase the noise level, thus leading 
to a decrease in the amplification since the dynamics of 
the particle is now hampered by the need to overcome a 
barrier. Therefore, the inflection point of the effective po- 
tential -which marks the appearance of the barrier- locates 
the position of a new optimal regime of noise in terms of 
signal amplification. 

2D numerical simulation. — In order to check the 
accuracy of the one-dimensional description we compared 
the obtained results with the results of Brownian dynamic 
simulations, performed by integrating the full overdamped 
Langevin equation (3). The simulations were carried out 
using the standard stochastic Euler-algorithm. 

The behavior of the spectral amplification 77 as a func- 
tion of the noise strength D for different tilting force values 
is depicted in fig. 5. It is worth to mention that the results 
of the ID modeling (lines) are in very good agreement with 
the numerical simulations of the full 2D-system (symbols) 
within a small relative error. As one would expect the res- 
onant behavior is absent for zero tilting force, i.e., for the 
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Fig. 5: (Color online). The spectral amplification 77 as a func- 
tion of the noise strength D at various tilting force values. The 
lines correspond to the numerical integration of the ID Fokker- 
Planck equation, eq. (17), whereas the symbols correspond to 
the results of the Langevin simulations for the two-dimensional 
structure with the shape defined by the dimensionless function 
w{x) = -0.25a;'* -HO. 5 3;^ + 0.01. 



purely symmetric case [20], whereas in the presence of a 
tilting force, the spectral amplification exhibits a double- 
peak structure. As discussed before, the inter-well dy- 
namics is responsible of the appearance of the main peak 
at higher noise strengths D whereas the second peak at 
smaller noise strengths is related to the disappearance of 
the barrier. By changing the tilt, the second peak gets 
shifted to higher values of noise following the predicted 
behavior for the inflection point, namely Dc ~ 0.144Fb. 
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Fig. 6: (Color online). Same as in fig. 5, but for various input 
signal frequencies. 

The behavior of the spectral amplification for different 
frequencies is plotted in fig. 6. The height of the main 
peak in 77 at high noise strength increases as the frequency 
of the input signal decreases, resembling the behavior of 
classical SR [2]. Overall, the perfect double-peak structure 
in the spectral amplification is present only at moderate 
frequency range. A similar behavior of double-peak SR 
one could also be observed in purely energetic systems, 
with a double- well potential, either at high input frequen- 
cies [27, 32] or at small frequencies in presence of inertia 
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Fig. 7: (Color online). Same as in fig. 5, but for various input 
signal amplitudes. 



effects [34]. 

Figure. 7 depicts the behavior of the spectral amplifi- 
cation for different values of the amplitude of the driving 
force. The inter-well dynamics responsible of the appear- 
ance of the main peak in 77 at higher noise strengths is not 
much affected by the variation of the input signal ampli- 
tude. In contrast, the second peak is greatly affected by 
the amplitude of the input signal. A clear double-peak 
structure in 77 can be obtained only at small signal ampli- 
tudes, since when the amplitude of the forcing exceed the 
bias Fb , the particle can move throughout the whole struc- 
ture assisted by the periodic forcing and therefore loses its 
sensitivity to the noise strength. 

Conclusions. — In this work we have studied the in- 
fluence of a constant bias on the resonant response of a 
geometrically constrained system. In absence of the tilt- 
ing force the system does not exhibit a SR-behavior [20]. 
Applying a constant bias, one can bring the system to 
an optimal regime where it exhibits a double ESR in the 
spectral amplification. We want to stress that the phe- 
nomenon is strictly entropic, since there is no purely ener- 
getic barrier in the system. The ESR emerges due to the 
interplay between the tilt and the entropic barrier which 
originates in our case from the unevenness of geometrical 
restrictions. The double-peak behavior arises due to two 
different dynamic regimes: (i) a regime, which is genuine 
of entropic systems, that occurs at high noise strengths 
where the main peak in the amplification is due to the 
synchronization of the periodic signal with noise-assisted 
hopping events; and (ii) a regime, at low noise strengths, 
at which the potential has an inflection point. The sys- 
tem thus becoming optimally sensitive to noise. Note that 
noise plays a two-fold role in an entropic system: on the 
one hand, it facilitates the dynamics, but on the other 
hand it intensifies the strength of the entropic contribution 
leading to the appearance of barriers that hamper trans- 
port. The second peak at smaller noise intensity signals 
the optimum balance between these two roles. This dou- 
ble resonant behavior renders the possibility for complete 
control of the optimization of the response of a system 
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in presence of noise. The main peak stemming from the 
inter-well dynamics can be strikingly controlled by tuning 
the frequency, while the second peak can be controlled ei- 
ther by changing the strength of the tilt or the amplitude 
of the periodic input signal. The double resonant peak 
behavior is a clear signature of the richness of dynamical 
behaviors in driven confined systems. The occurrence of 
this phenomenon depends only on a few controllable pa- 
rameters that can be advantageously tuned to achieve an 
optimal response. 
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